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The influence of a weak random potential on the collective modes of a trapped interacting Bose- 
Einstein condensate at zero temperature is calculated in the limit when the correlation length of the 
disorder is smaller than the healing length of the superfluid. The problem is solved in the Thomas- 
Fermi limit by generalizing the superfluid hydrodynamic equations to allow for the presence of weak 
disorder. We find that the disorder-induced frequency-shifts of the low-energy excitations can be of 
the same order of magnitude as the beyond mean-field corrections in the normal interaction recently 
observed experimentally. 



INTRODUCTION 

Some time ago, Huang and Meng [l[ have studied a 
homogeneous three-dimensional hard-sphere Bose gas in 
a random external potential as a model of superfluid he- 
lium in a disordered medium 0. In qualitative agree- 
ment with the experiments in porous media, they found 
that the formation of local condensates in the minima 
of the random potential reduces the superfluid compo- 
nent of the fluid even at zero temperature, where, in the 
absence of disorder, the whole fluid would be superfluid 
The recent experimental advances in trapping Bose- 
Einstein condensates in a disordered medium [j,|5| makes 
it interesting to test in a more quantitative manner the 
predictions of the model considered by Huang and Meng. 
For this purpose, we extend the latter approach to in- 
clude a harmonic trapping potential in addition to the 
weak external random potential. In the present work 
we consider a condensate in the limit of a large num- 
ber of particles N and in the presence of disorder with a 
correlation length shorter than the healing length of the 
superfluid. These conditions allow for a simple hydrody- 
namical formulation of the problem similar to the theory 
of wave propagation in random elastic media 0]. 

In the case of weak disorder, the corrections in the 
density profile or in the released energy of the Huang- 
Meng theory [l| turn out to be too small to be measur- 
able. Nevertheless, in the present paper we argue that 
disorder-induced shifts might be observable in the fre- 
quencies of the collective excitations 0, because these 
can be measured with an accuracy of a few fractions of 
a percent Q. More precisely, we show that long wave- 
length disorder can shift the frequencies of the low-energy 
excitations with the same order of magnitude but oppo- 
site sign as the beyond mean-field shifts due to repulsive 
atomic interactions Q in the gas. The latter effect has 
been recently observed experimentally by the Innsbruck 
group [lo| . Furthermore, in a harmonic trap, the influ- 
ences of disorder and interaction should be readily dis- 
tinguishable. This follows because, according to the gen- 



eralized Kohn's theorem [11, 12|, the interaction cannot 
alter the frequency of the lowest dipole mode 0] while, as 
we show, the latter is shifted by the presence of a weak 
random external potential. 

HYDRODYNAMIC EQUATIONS 

We consider a harmonically trapped Bose gas in an 
external random field. The grand-canonical Hamiltonian 
of the system is 



K- 



2m 



/i + [/ + Vc, 



(1) 

where ■(/'^(r), V'(r) are the field operators of an atom 
with mass m, fi is the chemical potential, and the in- 
teraction at low energy is described by the two-body 
T— matrix g = Airfi^a/m in terms of the s— wave scat- 
tering length a. The harmonic trapping potential is 
V^oxt (r) = m {ljjIx^ + ujyu'^ + bj'lz'^) /2. The oscillator 
frequencies define the harmonic oscillator length oho = 
(/i/mwHo)^^^, where who ~ {'^x'^y'^z f'^. The disor- 
der potential U is chosen with a Gaussian probabil- 
ity distribution characterized by the ensemble averages 
(C/(r)) = and {U{v)U{r')) = i?(r-r'), with R{v) = 
J c?k/(27r)'^e''"'i?(k). In the case of a fast-decaying 
disorder-correlation i?(r), the results of the theory do 
not depend significantly on its shape [l3|. We consider 
the case of a Gaussian correlation i?(k) = Roe~^ ^ 
where i?o and ^ characterize the strength and the corre- 
lation length of the disorder, respectively. 

A random potential in a Bose-Einstein condensate 
causes incoherent scattering which tends to localize and 
to deplete the condensate wavefunction. For a statisti- 
cally homogeneous system, Huang and Meng have shown 
[l[ that, even at T = 0, this results also in a depletion 
of the superfluid density Us . This latter must be distin- 
guished from the condensate density no and is related to 
the total density n by the relation Us = n — n„, where 
the normal (i.e. non-superfluid) component of the g 
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can be related to the disorder-induced depletion riR of 
the condensate through n„ = (4/3)nR [1}. Consistently 
with this picture, we assume that for weak disorder and 
at T = the superfluid component of the gas can be de- 
scribed by the phenomenological two-fluid hydrodynamic 
equations [1] 



d_ 

di 



n + V (VsTis -f v„n„) = 



"^^Vs + V (^1+ ]^m^fi ] = 0. 



(2) 



The hydrodynamical variables in this description are as- 
sumed to be ensemble averages over the realizations of 
the disorder potential. Therefore, the physical validity 
of Eqs. (21) must be considered as being restricted to the 
self- averaging regime where the wavelength of the hy- 
drodynamic modes is much larger than the length-scale ^ 
of the disorder potential. This condition determines the 
range of validity of our theory. The equations ([2]) can be 
understood as a Landau "two-fluid" model, for the un- 
pinned part of the condensate as superfluid component 
and the localized pinned condensate as normal compo- 
nent. We use these equations to describe the collective 
excitations of energy hw « fi^HO ^ A* in the presence 
of the harmonic trap. In the Thomas-Fermi regime and 
for long- wavelength oscillations IJ] , the non- uniform sys- 
tem can be considered locally as homogeneous, with the 
space-time-dependent density n (r, t) and the chemical 
potential related by /Li(r, t) =jM\n (r, t)] + V^^t (r) m the 
local density approximation [l4|. Here /i; [n (r, i)] is the 
chemical potential of a uniform gas at density n (r, t) . 
Because of the pinned character of the normal compo- 
nent, we must furthermore assume that only the super- 
fluid component reacts to the probe, while the pinned 
normal component remains stationary. This situation is 
reminiscent of the physically closely related problem of 
the fourth sound in ^He [1] and is expressed by the con- 
dition v„ = in Eqs. 0. 

Decomposing n (r, t) = n (r) 4- 5n (r, t) with the time- 
independent background density n (r) = Us (r) + n„ (r), 
and /X (r, = /ig + Sfi{r,t) with Sfi = (di-i/dn) 5n, and 
restricting ourselves to the linear regime, Eqs. ^ lead to 
the wave-equation 
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(r)V 



[n (r)] 
dn (r) 



Sn (r, t) 



=^ 0. 
(3) 



The disorder and the trap both appear implicitely in 
Eq. ^ via the superfluid density (r) and the equation 
of state fii [n (r)] which are determined in the following 
two sections, respectively. In the limit of a statistically 
uniform gas, Eq. ^ reduces to the wave-equation in a 
medium with a random refractive index 17 1. 



MEAN-FIELD THEORY IN THE TRAP 

In order to complete the equation of motion ([3]) we have 
to determine both the equilibrium superfluid density and 
the equation of state. To this purpose we now extend 
the Huang-Meng theory to a non-uniform system, by 
making use of the local density approximation. Expand- 
ing locally in plane waves \E'(r) = (l/^/V) X)k '^k(r)e*'"", 
the presence of the condensate is taken into account by 
setting ak(r) = {y/no{r)V)do^k + Sak{r). Retaining only 
the quadratic terms in the excitations 6ai^, fo^t from the 
condensate, the truncated Hamiltonian obtained from 
Eq. ([T]) can be diagonalized by the Bogoliubov trans- 
formation [H (5ak(r) = Uk (r) — t^k (r) a^j^ — Zk (r) 

and Sal^{r) = ""k l'^) '^k ^ '^k ('^) "^-k ^ ^wi^)' where 
the coherence factors Wk and iJk and the complex 
numbers Zk can be taken real by appropriately choos- 
ing the phase of the complex fields. Then we have 
ul (r) = {1 + [ek - /i + 14xt (r) + 2gn (r)]/Snk (r)}/2, 
fk (r) = '"k('') ^ 1' zi^{r) = 

[n (r) V]^ (7k K (r) - i^k (r)]' /hCii, (r), with 
Ek = h^k^/2m and, where the Bogoli- 
ubov spectrum [isj is given by ftfik (r) = 



Ek - M + Kxt (r) + 2.gn (r)] - [gn (r)] . The density 
is given by n (r) = fiQ (r) -I- riBog (r) -I- tt-r (r), where, be- 
sides the condensate density, (r) the local Bogoliubov 
depletion density [l^ JT-Bog (r) = (8/3)n (r) (r) o^/tt 
appears, and the depletion tir (r) induced by the random 
potential 



"R (r) = ^Ro (r) /R[47rn (r) ^^a] 



(4) 



Here, riR^ (r) = i?o (m^/S tt^''^ h'^) yjn (r) /a denotes the 
local depletion in the limit of a 5-correlated disorder po- 
tential and the function lldl 



Ux) 



(1 + 4a;) Erfc [^2^^ - 2^2x/tj: (5) 



(with /r,(0) = 1) includes the effects of its non-vanishing 
correlation length ^. With tir (r) now in hand, the local 
superfluid density ns(r) is given in terms of n(r) by the 
relations ^^(r) n(r) — ?T.„(r), where the normal compo- 
nent is n„(r) = (4/3)nR(r) as discussed before. 

The total disorder-induced depletion of the conden- 
sate A^R = / drriR (r) , can be calculated to leading 
order in i?o by replacing in Eqs. (g]) the density n(r) 
by its zero-order Thomas-Fermi approximation n (r) ~ 
[mtf - K=,t (r)] /5, where //tf = ^{^f^ is the 
mean-field chemical potential of the Gross-Pitaevskii 
theory. In the case of (S-correlated disorder, integrat- 
ing over the mean-field Thomas-Fermi radius we get 
N^JN ~ (157r/32)[7iR„ (0)/n(0)]. The theory is valid 
when Njig/N <C 1, which can also be rewritten as 
R'q (r) = m-^i?o/87r^/2/i*A/n(r)a ^ 1, or, in local form 
as the condition riR^ (r) ^ n (r) . In the Thomas-Fermi 
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regime the latter condition is satisfied everywhere in the 
trap with the exception of a narrow shell at the bound- 
aries where the condensate density vanishes. In the 
present work we neglect possible small corrections from 
this effect. 



BEYOND MEAN-FIELD EQUATION OF STATE 

The ground-state density n (r) at equilibrium in 
Eq. Q can be calculated by using the local density ap- 
proximation for the chemical potential /ip = Mi ['^ (r)] + 
Vcxt (r)- In the mean-field Huang-Meng theory outlined 
in the previous section, the latter is determined by the 
Thomas- Fermi result //tf = 5'^(r) + 14xt (r) of the Gross- 
Pitaevskii theory. In this mean-field app roximation the 
effects of quasi-particle interactions [20| as well as the 
scattering between the quasi-particles and the impurities 
[lit neglected. In order to incorporate these processes 
in the calculation of the collective modes, the ground- 
state value of the chemical potential fi has to be fixed in 
terms of N — J drn(r) by including beyond mean-field 
corrections in the equation of state. These latter can be 
calculated in the framework of the Bogoliubov theory [l3| 
resulting in 



m (n) ^ng 




(6) 



where /iR,; (n) = Ggn^^ (r) /^[47rn (r) ^^^j^ f^(^^^ 



e^^ (3 + Ax) Erfc (V^x) - 2^/2^ /3. In the following 
we wish to focus on the effects of disorder. Therefore, 
in the equation of state ([6]) we can neglect the beyond 
mean-field corrections due to the normal interactions, 
keeping only those due to the disorder. This is possible 
because, although the two different corrections can be of 
the same order, they lead to effects which are additive to 
lowest order. Moreover, in an experiment one could focus 
on the effects of the disorder by tuning, via a Fcshbach 
resonance, to the regime R'q{0) ^ y^n{0)a^ where they 
become dominant. Other corrections due to finite-size, 
non-linearity and temperature have been sufficiently dis- 
cussed in Ref. in connection with the frequency shifts 
induced by the be yon d mean-field effects in the theory of 
Huang and Yang [20|. 

Using Eq. ([6|), for small disorder, we can find the equa- 
tion for the ground-state density by iteration. Defining 
the Thomas-Fermi density as utf (r) = [a'o ^ Kxt (r)] /g, 
with Ho{N) determined from the normalization condition 
N = J drn{r) including the disorder correction in the 
equation of state ([6]), we find 

n (r) ~ riTF (r) - 6 n^, (r) ^ [47rnTF (r) ^'a] . (7) 

In the derivation we have simplified the disorder-induced 
depletion tirq (r) of Eq. ^ via replacing n by its zero- 
order approximation tt-tf- The same simplification is 



made in the argument of /r and will be used in that 
of the function defined below. 

Another quantity we need to evaluate in the hydrody- 
namic equation ([3]) is the term proportional to the inverse 
of the compressibility. From Eqs. ([U and ^ we have 



(dfii/dn) =g{l + 3[nn„ (r)/nTF (r)]/^ [47rnTF (r) C'a] 
+ 6[nR„ (r) /nxF (r)]/^ [47rnTF (r) ^^a] }, 



(8) 



with/x(x)=x |e2=^ (5 + 4x) Erfc (\/2^) 



4 / 2 1+x 

3 Y TT ^/k 



COLLECTIVE MODES 

Having calculated the beyond mean-field corrections 
due to disorder in the equation of state, we can proceed 
to determine the explicit form of the hydrodynamic equa- 
tion describing the low-energy collective modes of the 
system in the linear regime. Using the result of Eq. ([7]) 
and Eq. ([8]), and retaining only terms linear in i?o, the 
hydrodynamic equation of Eq. ([3]) can be put into the 
final form 



m-^Sn—\/ [guTp'^ Sn] = — V 



4nRo/R 



g\76n 



(9) 



{g [3 nRo/^ - Qnnofxl Sn} , 



where we have left the r— dependence of the coeffi- 
cients implicit. In the limit of a uniform gas the so- 
lutions of Eq. ([9]) exhibit a phonon dispersion fiuj = 
cq. Using Eq. ^ and expanding for ^ <ti Chcai we 



find 



which describes the 



shift of the velocity of sound induced by the disorder 



from the Bogoliubov result Cq = 
system with healing length ^hcai 



gn/m of the clean 
= l/VS-Trna. When 
[13, El 



^ = this reproduces the result found in Refs 
for 5— correlated disorder. Moreover, putting the r.h.s. 
of Eq. ([9]) equal to zero we recognize the mean-field 
equation used by Stringari in Ref. [3| in the case of 
an isotropic trap. This yields the dispersion relation 

1/2 

cJo(?^r,0 = ujYio i2nl + 2nrl + ?>nr + t) , for excita- 
tions with rir radial nodes and I angular momentum. In 
the presence of disorder, Eq. ([9]) can be solved consider- 
ing the right-hand side as a small perturbation. Defin- 
ing the function h[r) = 3 riR^, (r) [47r^^nTF (r) «] — 
6??R(, (r) [47r5^nTF (r) a] , we find ultimately the fre- 
quency shift 

'^^o(e)_ .9 S dY[{V^5n*)h5n + 5n*V {^n^J^ySn)] 



2muJn 



J dr 5n*5n 



(10) 



where dn are the solutions of the mean-field equation. In 
contrast with the beyond mean-field effects due to the 
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atomic interaction described in Ref. the frequencies 
of the so-called surface-modes are influenced by disorder. 
These modes have the principal quantum number n^. ~ 
and represent solutions of the type (5n (r) ~ r^Yim with 



a mean- field dispersion law given by [7| cJo('^r = 0) = 
^HoV^- Because they satisfy the condition V^Jn (r) = 
the left term in the numerator of the r.h.s. of Eq. pO|) is 
zero. Note that this term comes from the change in the 
macroscopic compressibility contained in the term dfi/dn 
in the equation of state of Eq. ([3|). However, the right 
term carrying the factor 4/3 does not vanish. This cor- 
rection originates because the normal part of the liquid 
remains stationary under an external probe. In general 
we obtain for the frequency shift of the surface modes 



5uo{£.) _2l + 3 
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R'o{0)Io,l,n{i) (11) 



with 



lo.Lni I n = 2 



or 



independent from the quantum number m. The di- 
mcnsionlcss variable ^ is given by the relation = 

($i?TF/aHo) /2- III the case of the dipole mode 5n ~ 
rcoa9 with I — 1, m = 0, the mean-field result 
Wdip = i^HO coincides with the harmonic oscillator result 
ciJosc = t^HO (2?T,r -l- 0- This follows from the fact that 
in a harmonic potential the lowest dipole mode (n,- — 0, 
/ = 1) corresponds to the oscillation of the center of mass, 
and is unaffected by the interatomic forces (Kohn's the- 



orem) However, during its motion, the supcrfluid 

density is "hampered" by the normal component which 
remains stationary. In the limit of (5— correlated disor- 
der, Eq. In]) gives Siudipi^ = 0)/wdip = -(57r/16)i?(, (0). 
In the same limit, the shift of the quadrupole mode 
I = 2, m = 2, described by 6n ^ rsin^fle^""^, is 
6luq{^^0)/ujq = -(357r/96)i?() (0). According to the 
definition of the dimensionless parameter R'q (0), we see 
that, besides disorder, also the presence of interactions 
and thus of superfluidity, is crucial in order to have the 
effect. 

The effects of a non-zero disorder correlation length ^ 
on the frequency shift of the surface modes can be better 
understood when considering the ratio 

Scja{£_) 2 r{l + 2) 



6uJo{0) y/^T{l + 3/2) 



(12) 



obtained from Eq. pip . In Fig. [T]the result is illustrated 
in the case of the dipole and the quadrupole oscillations. 
We plot the relative frequency shift Sujo{0 / ^^o{0) as a 
function of the dimensionless variable = ■C^tf/QhoV^ 
in the case of the surface excitations with I = 1 and 
1 = 2. The shift decays rapidly with increasing ^. How- 
ever, as we have anticipated above, the approximation 




Figure 1: Relative frequency shift 5uJo{i,) / 5uJo{0) as a function 
of the variable ^ = ^Rtf / tSHO in the case of the surface 
excitations with I = 1 (solid line) and I = 2 (dashed line). 



introduced in order to derive the hydrodynamic equa- 
tion limits the range of validity of our theory to 
random potentials with coherence length ^ much smaller 
than the heahng length ^heai = 1/ \/8Trn{0)a of the su- 
pcrfluid. By using the relation between the central den- 
sity n(0) and the Thomas-Fermi radius we have that 
6ioai = -RTF(aHo/-RTF)^ arid the inequality ^ < ^hoai 
can thus be rewritten as ^ 1- Therefore, the results 
of Fig. [1] for the region V2C ^ 1 can only be considered 
as an extrapolation. 

Experimentally, the effects of a random potential on 
the dipole and quadrupole modes of a trapped Bose- 
Einstein gas have been investigated in Ref. [4| by using 
optical laser speckles as realizations of disorder configu- 
rations. In these experiments, the smallest length scale 
of the speckle potential is of about ^ ~ 10 /xm while the 
axial Thomas-Fermi radius is 100 /xm. As a typical ex- 
perimental situation where the Thomas-Fermi condition 
Na/ano 3> 1 is satisfied, we can assume i^xF/oHO — 10- 
This would imply a value ^ ~ 6 in our description, which 
is far beyond its range of validity. In particular, for 
large values of the disorder correlation length, such that 
C ^ Cheai, the system is not self- averaging and the ob- 
served frequency-shifts must depend on each individual 
realization of the speckle potential [i^, [l^- In that case, 
the disorder average must experimentally be taken by 
determining, for each mode, the mean value of the mea- 
sured frequencies over many different realizations. Both 
experiment Q and theory predict no shift in that 
regime in qualitative agreement with our extrapolation. 
Decreasing further the correlation length of the disorder 
would allow to enter the regime f ^ Cheai, where the 
system becomes self-averaging, and the shift predicted in 
our theory should become observable. 

Differently than for the surface modes, the mean- 
field frequency ojm — \/5who of the lowest compres- 
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sional mode (n^ = 1, / = 0) is shifted by the beyond 
mean-field corrections in the normal interaction. The 
shift has been calculated in Ref. [§] and is SwmI^m = 
(637r/128) ^/a?n{^. Inserting the density oscillation 
6n (xl — (5/3) the analogous correction due to the dis- 
order calculated from Eq. pH]) in the limit = amounts 
to (5wm/wm = -(4697r/768)i?;, (0). 



ANISOTROPIC TRAP 

The results can be generalized to non-sphcrical traps, 
by considering an axially deformed trap of the form 
Kxt (r) = m [uj\r']_+ uj^z"^) /2, where r_L = \J + 
is the radial coordinate. In order to put better in ev- 
idence the role of the anisotropy, we restrict ourselves 
to the case of disorder with vanishing correlation length. 
In the presence of anisotropy, the functions of the form 
8n ^ r^Yim are still solutions of the l.h.s. of Eq. ([5]). In 

1, m = 0) with 
tJz, the disorder 



the case of the dipolc oscillation (/ 
the mean-field mode frequency o^dip 
corrections lead to the same relative shift Jwdip/i^dip as 
in the case of the isotropic trap. The same argument 
applies to the quadrupole mode / = 2, m = 2, with 
= v^^x- The quadrupole mode with m = involves 
a mixing with the monopole mode = 1, i = 0. In this 
case the dispersion law is [wo^(™ = 0)]^ = ^1(2 + 
fA^ ± iVSA-* - 16A2 + 16), where A = wjw^ character- 
izes the deformation of the trap. Using the corresponding 
oscillation 8n ^ -^{{^f - 2] + ri + \{^f - 4].^ 
in Eq. (dH]) we find the shifts 



-7^^o(0) 



±72 ± 9A2 -f 107\/T6~T6A2+W 



1536V16 - 16A2 + gA-i 



(13) 



Note that for a spherical trap (A — 1) we have 
lOq = -s/Swho and lOq = v^who and Eq. (fT3|) re- 
covers the shift of the quadrupole mode JtJq/tjQ = 
— (357r/96)i?Q (0) and the shift of the monopole mode 
iJtJM/wM = — (4697r/768)i?o (0) calculated above. The 
A— dependence of the frequency shifts of Eq. (fT^ is shown 
in Fig. [5] by plotting the functions (5ijj^(A)/(5ijJm (lower 
curve) and Jwj^ (A)/(5a)Q (upper curve). 



CONCLUSIONS AND OUTLOOK 

We have calculated the shifts to the collective frequen- 
cies of a zero-temperature trapped Bose gas induced by a 
weak external random potential with a correlation length 
smaller than the healingh length of the superfluid. We 
have shown that the realization of such a limit in cur- 
rent experiments in trapped Bose-Einstein condensates 
in random media would allow to measure for the first 
time the predictions of the Huang and Meng theory in 




Figure 2: Relative frequency shift (Jo^q" (A)/(5[^m (lower curve) 
and iijjQ (A) /iSctJQ (upper curve) of the m = modes, result- 
ing from the coupling between the quadrupole and monopole 
modes, as a function of the deformation parameter A = 
LJz/uj±. The disorder correlation length is taken here C = 0. 



a quantitative way. Moreover, the interplay between the 
pinned and the unpinned components of the condensate 
is expected to produce a deviation from the generalized 
Kohn's theorem for the center of mass motion of the su- 
perfluid. This latter phenomenon provides an unambigu- 
ous signature of the disorder-induced effects against the 
beyond mean-field corrections due to the interatomic in- 
teractions. 

These results have been derived by means of a hy- 
drodynamic approach and could be extended to consider 
the strongly anisotropic traps which realize the so-called 
"one-dimensional mean-field" regime described in Ref. 
24j . This would be relevant in connection with the recent 
experimental and theoretical investigations of the trans- 
port of Bosc-Einstcin condensates in one-dimensional mi- 
crotraps in the presence of disorder [2l,[2i,[23,[2i,[2^. 
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